Comparison of Dynamical Heterogeneity in Hard-Sphere and Attractive Glass 
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Using molecular dynamics simulations, we have determined that the nature of dynamical hetero- 
geneity in jammed liquids is very sensitive to short-ranged attractions. Weakly attractive systems 
differ little from dense hard-sphere and Lennard- Jones fluids: Particle motion is punctuated and 
tends to proceed in steps of roughly a single particle diameter. Both of these basic features change 
in the presence of appreciable attractions. Transient periods of particle mobility and immobility 
cannot be discerned at intermediate attraction strength, for which structural relaxation is greatly 
enhanced. Strong attractions, known to dramatically inhibit relaxation, restore bimodality of par- 
ticle motion. But in this regime, transiently mobile particles move in steps that are significantly 
more biased toward large displacements than in the case of weak attractions. This modified feature 
of dynamic heterogeneity, which cannot be captured by conventional mode coupling theory, verifies 
recent predictions from a model of spatially correlated facilitating defects. 



Dynamical heterogeneity is perhaps the most reveal- 
ing feature of relaxation in deeply supercooled liquids. 
A high-temperature fluid is dynamically homogeneous in 
the sense that the local environment restricting fluctua- 
tions of any given particle is, for all important purposes, 
identical to that surrounding any other particle, even 
on the short time scales of basic microscopic motions. 
The distribution P(8r, t) of particle displacements <5r as a 
function of time t provides a quantitative measure of such 
uniformity. Results of molecular dynamics simulations 
indicate that P(6r, t) is Gaussian over a wide range of dis- 
placements, as would be expected from a mean-field per- 
spective, for typical dense fluids [1|. The microscopic en- 
vironments constraining particle motion in a glassy ma- 
terial are by contrast profoundly nonuniform, even on the 
long time scales of large- wavelength relaxation[2],|3j,U,|5|. 
This fact has been clearly demonstrated by experiments 
that focus on dynamics of single probe molecules |6J or 
subsets of molecules in a pure liquid that relax more 
slowly than the average [jj |S|- As a result, P(Sr,t) de- 
velops substantial weight in the wings, reflected in ap- 
preciably nonzero values of the non-Gaussian parameter 
a 2 (t) = |((5r 4 (i)}/((5r 2 (i)} 2 - 1. Microscopy studies of 
colloidal suspensions have confirmed this expectation|9|. 

Through extensive computer simulations, a detailed 
picture of dynamical heterogeneity in simple jammed liq- 
uids (e.g., a binary mixture of Lcnnard-Jones spheres) 
has developed[lfJ. At low temperatures, the majority of 
particles are confined within cages, composed of neigh- 
boring particles, that may persist for very long times. 



Eventually, a rare, collective rearrangement frees a parti- 
cle from its cage, transiently allowing it to move rapidly 
over distances comparable to a particle diameter. Since 
such a displacement itself facilitates local rearrangement, 
transient mobility appears to propagate continuously and 
with some degree of directionality. Schematic models of 
glassiness have been constructed with only these features 
in mind 11, 12, 13]. They account for a surprising va- 
riety of anomalous behaviors and yield unique scaling 
predictions for the length and time scales characterizing 
relaxation|13l Il4i llo| . 

These basic features of dynamical heterogeneity have 
little to do with the identity of particles comprising a 
glassy material or with the interactions between them. 
They have been reported for many model atomic liquids 
as well as for viscous silica, a network-forming liquid 
which vitrifies with qualitatively different temperature 
dependence 16]. It is therefore tempting to presume that 
the scenario sketched above is universal among super- 
cooled molecular liauids [12lll7| . 

When interactions between particles in a simple liquid 
are augmented by strong, short-ranged attractions, spa- 
tially ave rag ed dynamical quantities change in nontrivial 
ways [la llil l2ii l21| . This situation has been realized 
experimentally by adding linear chain molecules to sus- 
pensions of colloidal particles, whose direct interactions 
are almost purely repulsive [22|. By varying <f> p one can 
tune the fluid from purely repulsive to strongly attrac- 
tive. For large colloid volume fraction, </> c , and vanishing 
attraction strength (<f> p = 0), the suspension is in essence 



a dense hard-sphere fluid with the basic phenomenology 
of simple supercooled liquids. As <p p increases, however, 
relaxation accelerates significantly, so that a hard-sphere 
glass can be "melted" by adding attractions 18] . Be- 
yond a certain value of (f> p , relaxation becomes instead 
more sluggish as more polymer is added, leading to re- 
vitrification at large <fi p . 

In this paper we examine whether the peculiar behav- 
ior of attractive colloids reflects fundamental changes in 
the nature of dynamical heterogeneity. Although exten- 
sive computational work has been done to characterize 
such heterogeneity in simple liquids, previous simulations 
of materials with short-ranged attractions have not fo- 
cused on correlated microscopic motions. For this pur- 
pose we have adopted the model of Puertas et al. for 
polymer-mediated interactions between colloids |l3, |2fJ. 
The effective interaction potential between a pair of col- 
loids separated by distance r, described in detail in Refs. 
21-22, is plotted in Fig. ^ for the values of <p p we have 
studied. Interactions consist of a short-ranged repulsion 
parameterized by the sum of particle radii, a short-ranged 
attraction that mimics the polymer-induced depletion in- 
teraction, and a slowly varying, long-ranged repulsion 
designed to prohibit phase separation [23. The effective 
range and form of each term are precisely as given in Refs. 
21-22. We have focused on weakly polydisperse systems, 
in which these radii are drawn from a uniform distribu- 
tion, p(R) = (2 5a)~ 1 9[R-(a-5a)]9[(a + Sa)-R], with 
mean a and half-width 5a = a/10. The Hcavisidc funcion 
8(x) is defined as usual to be 1 for x > 0, and for x < 0. 
In this paper all quantities with units of length have been 
scaled by the mean radius a, and quantities with units 

of time have been scaled by t = yfy-- We have used 
standard methods of molecular dynamics to propagate 
equilibrated systems of 1000 periodically replicated col- 
loidal particles, stochastically rescaling colloid momenta 
every 101 time steps to maintain a Boltzmann distribu- 
tion at temperature T = | . We have investigated a single 
colloid volume fraction </> c = ^0-a 3 (l + (Sa/a) 2 ) = 0.55 
and several polymer volume fractions (i.e., attraction 
strengths as plotted in Fig. QJ ranging from (f> p = 0.05 
to 4> p = 0.375. 

In order to contrast features of dynamic heterogeneity 
unique to attractive systems with those generic to sim- 
ple supercooled liquids, we have selected a colloid den- 
sity for which relaxation is already very slow for 4> p = 0. 
The self-diffusion constant D((f> p ), plotted in Fig. Q] as a 
function of (j) p , thus evinces the re-entrant behavior we 
have described. Specifically, D increases by two orders 
of magnitude as <f> p approaches 0.25 and then decreases 
sharply for higher polymer concentrations. Any mean- 
ingful comparison of relaxation mechanisms at different 
cj) p must take into account the disparate time scales cor- 
responding to this range of diffusivity. Here we examine 
time evolution over intervals scaled such that the overall 
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FIG. 1: (a) Diffusion constant D as a function of polymer 
volume fraction <j> p . In each case, D was computed from the 
mean squared displacement averaged over particle size. In- 
set: Potential energy of interaction V(r) as a function of the 
distance r between two particles with radius a for the values 
of (j>p we have studied numerically. From the least negative 
minimum of V(r) to the most negative minimum, these val- 
ues are <f} p = 0.05, 0.15, 0.2, 0.25, 0.3, 0.35, and 0.375. (b) 
Non-Gaussian parameter 0:2 (£) as a function of time t. 



extent of relaxation is comparable for each state. We use 
several different measures of the extent of relaxation, in- 
cluding a%{t), mean squared particle displacement, and 
the dynamic structure factor. Although these choices 
are somewhat arbitrary, they paint a consistent picture 
of changes in dynamic heterogeneity induced by short- 
ranged attractions. 

The non-Gaussian parameter, plotted in Fig. \$£b) as 
a function of time for several values of </> p , indicates 
that short-ranged attractions effect changes more pro- 
found than simply a renormalized average time or length 
scale |l9ll20| . The peak of oti(t) roughly locates the time 
of maximum dynamic heterogeneity. We denote the time 
corresponding to this peak as t* . The dependence of t* 
on (ftp closely mirrors that of the diffusion constant, re- 
flecting a global change in relaxation time. But the shape 
and scale of a-i (t) also change significantly with attrac- 
tion strength. Most notably, the peak height a.2{t*) de- 



clines by nearly an order of magnitude as <j) p approaches 
0.25, then grows rapidly for larger values of (j> v . This 
evolution strongly su gges ts a change in the character of 
microscopic dynamics [l9j. 

The full distribution of particle displacements, or of 
their logarithm (P[log 10 (|<5r|), £]), over a specific time in- 
terval (of duration t) provides a more detailed picture 
of microscopic rearrangements J24| . Fig. [5] shows a plot 
of P[log 10 (|<5r|),To] for states representative of weak at- 
tractions (0 P = 0.05), strong attractions (<j) p = 0.375), 
and intermediate attraction strength (</> p = 0.25). Since 
relaxation rates vary greatly a mong these three states, 
we have followed Cates et al. ,24) in comparing mo- 
tion over time intervals tq(4> p ) yielding the same mean 
squared displacement, 10a 2 . As reflected by a2(t), the 
distributions at cf> p — 0.05 and (j> p — 0.375 are highly 
non-Gaussian, exhibiting distinct populations of espe- 
cially mobile and especially immobile particles. Cates et 
al. |24j have reported a similarly bimodal distribution of 
log(|<5r| 2 ) for strongly attractive colloids at much lower 
densities (cj) c = 0.4) and have suggested that such bi- 
modality is a unique feature of attractive glasses. Multi- 
peaked van Hove distribution functions, however, have 
been reported several times for systems lacking short- 
ranged attractions |25l 12(1 |27| . Indeed, Fig. 2 demon- 
strates that the distinction between mobile and immobile 
particles is in fact more pronounced when attractions are 
almost negligibly weak. Bimodally distributed particle 
displacements thus appear to be a feature common to 
many sluggish systems. 

The glassy states with weak and strong attractions also 
share a degree of structure in P[log 10 (|<5r|), t ] within the 
subpopulation of mobile particles. For 4> p = 0.05 parti- 
cles clearly tend to move in discrete steps of approxi- 
mately integer multiples of a typical particle diameter, 
2a. This feature highlights the decoupling of diffusion 
and structural relaxation in jammed liquids. Although 
fluctuations in local environment may permit a particle 
to move out of its cage, density correlations persist such 
that the newly formed cage has a well-defined spatial re- 
lationship with the original. Stepwise motion is much 
less pronounced at high (f> p . 

By contrast, the shape of P[log 10 |(5r|,r ] at <f) p = 0.25 
is nearly that of a Gaussian, plotted for reference in 
Fig. [21 There is evidence neither of a distinct popula- 
tion of immobile particles nor of stepwise motion. The 
statistics of single-particle displacements at intermediate 
attraction strength thus more strongly resemble those of 
a dynamically homogeneous fluid at lower density than 
those of a non-attractive fluid at the same density. The 
simplest conclusion is that the <j) p = 0.25 system is 
in essence dynamically uniform. Evidence exists, how- 
ever, that correlated motions of neighboring particles do 
show signs of dynamic heterogeneity 28] . This situation 
could be expected if particle displacements were domi- 
nated by movement of clusters transiently stabilized by 




FIG. 2: Distributions P[log 10 |<5r|, To] of the logarithm of par- 
ticle displacements |<5r| at time to. Results are shown for 
polymer volume fractions <j) p = 0.05, 0.375, and 0.25, and for 
a Gaussian distribution of or. Times To((j>) and the width of 
the reference distribution were chosen to obtain a consistent 
mean squared displacement, (|o~r| 2 ) = 10a 2 . 



attractions 28J . Heterogeneity associated with formation 
and decay of a cluster would not have a strong signature 
in P[log 10 |<5r|,To] due to translation of the cluster as a 
whole. 

The essence of a dynamic heterogeneity perspective on 
glassy liquids is that relaxation over any short interval is 
driven by a small subset of particles that are temporarily 
much more mobile than the average. The statistics of 
extreme displacements should therefore be revealing of 
basic relaxation mechanisms 10]. Here we focus on parti- 
cles among the 5% most mobile over an interval of length 
t* C To. We judge mobility in this case by monitoring 
the largest displacement |<5r| max of a particle from its po- 
sition at the beginning of each interval. Distributions of 
maximum displacement magnitudes for these especially 
mobile particles, P^Hl^lmaxji*); are plotted in Fig. 
for each of the polymer volume fractions we have stud- 
ied. For purposes of comparison, we have scaled |<5r| max 
by its most likely value ro((f> p ) for each <f) p . 

For reference we have included in Fig. [3] an extreme 
value distribution 7-g ( | <5r | ) that would be obtained for 
a Brownian analogue of our system. Because maximum 
displacements are not easily defined for fractal trajecto- 
ries, we consider in this case the displacement of a par- 
ticle \5t\ from its initial position only at the end of an 
interval. The interval length is in fact arbitrary, deter- 
mining only the overall scale of displacements, which is 
irrelevant for the comparison in Fig. [31 In detail, we com- 
puted Pg (|<5 r l) by repeatedly drawing N = 1000 dis- 
placements from a Gaussian distribution with unit vari- 
ance in all three dimensions, each time adding the 0.05-/V 
largest values to a histogram. Pg (|<^ r l) is therefore a 



superposition of extreme value distributions: 
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Here, <7j(|r5r|) is the probability density for observing \Sr\ 
as the j th largest value in a sample of size N. In the limit 
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where Xj — a/2 In (j /N). Because convergence to this 
asymptotic limit is slow, however, we have chosen to con- 
struct Vft (|<5r|) from direct sampling. 

The scaled distributions in Fig. [3] reveal a monotonic 
trend toward broadly distributed mobile particle dis- 
placements as attraction strength increases. For weak at- 
tractions, the large-displacement tail of P^' (\Sr\ max , t*) 
is attenuated relative to a simple Brownian fluid. For 
strong attractions, this tail is greatly enhanced. Al- 
though the overall shape of the full displacement dis- 
tribution for intermediate attraction strength is roughly 
Gaussian, statistics of the extreme subensemble demon- 
strate that mobile particles nonetheless execute larger 
jumps (relative to the average) than in a Brownian refer- 
ence system. This fact is consistent with the attraction- 
enhanced large-displacement tail of P[log 10 |<5r|, To] plot- 
ted in Fig. El 

The comparison in Fig. [21 emphasizes changes in the 
shape of displacement distributions. Because we have 
scaled distance by different lengths for different values 
of <j)p , these results do not directly imply that transiently 
mobile particles execute larger jumps in space in the case 
of strong attractions. They instead reflect the breadth 
of particle motion relative to the mean, a feature that 
could reflect tightening of particle cages as well as growth 
of facilitated regions. More direct evidence for both of 
these trends is provided by scrutinizing the trajectories 
of individual particles. 

The routes traced by several particles are plotted in 
Fig. 01 for the smallest and largest values of p . We have 
chosen the duration of displayed trajectories to be r Q , 
the time required for correlation of density fluctuations 
over length scales comparable to a particle diameter to 
decay by a factor of e. Dynamics over this time scale, 
which is much larger than t* in both cases, exhibit the 
full character of particle motion. This choice also al- 
lows a straightforward comparison of dynamics for weak 
and strong attractions, since the corresponding values of 
r a are very similar (r Q w 8 x 10 3 for </> p = 0.05 and 
r a « 10 4 for 4> p — 0.375). The depicted excursions of 
extremely immobile (top panels of Fig.Q and extremely 
mobile (bottom panels) particles reinforce the dynamical 
features we have gleaned from probability distributions. 
The 5% least mobile particles exhibit only small flucuta- 
tions about their initial positions, both for p = 0.375 (a) 
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FIG. 3: Distributions -P^'d^rlmax, £*) of maximum displace- 
ments for the 5% most mobile particles. We use as a measure 
of mobility |r5r| max , the largest displacement a particle under- 
goes (from its position at time t) during the time interval t to 
t + 1* . Distributions and displacements have been scaled to 
coincide at the peak of P^ y '. A scaled extreme displacement 
distribution is also plotted for a perfectly Brownian system of 
1000 particles. 



and for (f> p = 0.05 (c), even over this long time scale. The 
cages that constrain these particles' motion are clearly 
smaller in the case of strong attractions. The 5% most 
mobile particles, on the other hand, move several parti- 
cle diameters during the same interval. For the case of 
weak attractions in panel (d), the discrete nature of par- 
ticle motion and correlations between subsequent cages 
are immediately evident. Trajectories are qualitatively 
different for the case of strong attractions in panel (b). 
Here, mobile particles explore much more diffuse regions. 
Domains of facile particle motion are clustered in space 
and markedly elongated, with asymmetries apparently 
correlated over several particle radii. 

Most of the qualitative changes in dynamical hetero- 
geneity we have reported can be understood as con- 
sequences of changing spatial patterns of structural 
defects 31]. Models based on dynamical heterogeneity 
typically assume that the subtle defects which enable lo- 
cal relaxation are sufficiently sparse as to be statistically 
independent (despite significant correlations in defect dy- 
namics). We have proposed that an important effect of 
short-ranged attractions is to introduce non-negligible 
spatial correlations among such facilitating entities[30J. 
In particular, defects should aggregate with increasing 
attraction strength as an indirect result of particle clus- 
tering. Microscopic regions of mobility thus grow in size, 
but become still more sparse if their overall concentra- 
tion (loosely analogous to free volume) is held fixed. This 
picture accounts for the broadening distributions of mo- 
bile particle displacements we have computed. Particle 
trajectories depicted in Fig. 0] make the agreement es- 
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FIG. 4: Representative trajectories of particles whose dis- 
placements in the xy-p\ane over a time interval r a are much 
smaller or much larger in magnitude than average. Particle 
traces have been projected onto this plane for graphical sim- 
plicity. Left panels depict dynamics of a strongly attractive 
system, <j> p = 0.375. Particles among the 5% least mobile are 
shown in (a), while particles among the 5% most mobile are 
shown in (b). Right panels correspond to a weakly attractive 



system, 



0.05. Again, trajectories exhibiting extremely 



small displacements (c) are plotted in the top panel. Trajecto- 
ries exhibiting extremely large displacements (d) are plotted 
in the bottom panel. 



pecially vivid. The limited spatial extent of facilitating 
defects in a hard sphere glass cuts off the range of avail- 
able displacements. Clustering of these defects as at- 
tractions are introduced provides increasingly extended 
loose regions for mobile particles to explore. Our nu- 
merical results provide strong evidence for segregation of 
jammed and unjammed regions of attractive liquids at 
high density. Such segregation is an obvious feature of 
attractive colloids at low <j) c , which form stable (though 
likely noncquilibrium) gel-like networks. In that case, the 
spacing between dense regions of the network establishes 
a minimum length scale for dynamic heterogeneity. It 
is remarkable that remnants of this behavior persist at 
high packing fraction, where spatial heterogeneity of liq- 
uid structure is subtle. 

In summary, the change in dynamics induced by short- 
ranged attractions in a dense model liquid is dramatic 
even on the microscopic scale. Our results reveal three 
distinct regimes of dynamical heterogeneity. For weak 
attractions, mobilized particles make discrete jumps be- 
tween cage structures, which may remain correlated over 
many jump times. For a range of intermediate attrac- 
tion strengths, Gaussian particle displacement statistics 



suggest instead very fluid and uniform motion. The role 
of attractions in this regime, we suggest, is to bind small 
transient clusters which move on a time scale comparable 
to their lifetimes. Strong attractions restore some dis- 
creteness of particle displacements, presumably because 
transient clusters are too large to move appreciably on 
pertinent time scales. We thus conclude that "attractive" 
glassiness is driven by spatial redistribution of facilitating 
defects. Although extended loose domains permit large 
particle displacements, their growth depletes mobility in 
surrounding areas, which in turn inhibits relaxation of 
domain interfaces. We are pursuing further calculations 
to confirm this clustering of mobility in dense environ- 
ments and to characterize the correlated fluctuations un- 
derlying fluidity at intermediate attraction strength. 
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